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l/~) ' Subject of this work are the applications of a field theoretical model, called here 

generalized nonlinear sigma model or simply GNLerM, to the dynamics of a chain 
subjected to constraints. Chains with similar properties and constraints have been 
discussed in a seminal paper of Edwards and Goodyear using an approach based on 

a ■ 

the Langevin equation. 

The GNLctM has been proposed in a previous publication in order to describe 

a: 

the dynamics of a two dimensional chain. In this paper the model is extended to 

l o j . 

d dimensions and a bending energy term is added to its action. As an application, 

(N ■ 
£> ■ two observables are computed in the case of a very stiff chain. The first observable 

*n : 

is the dynamical form factor of a ring shaped chain. The second observable is a 

°' 

straightforward generalization to dynamics of the static form factor. This observable 

t> ' 

^) • is relevant in order to estimate the average distance between two arbitrary points of 

oo : 

<0 ■ the chain. 

.£h \ Finally, a variant of the GNLcM is presented, in which the topological conditions 

>< ■ 

which constrain the motion of two linked chains are imposed with the help of the 
d | 

Gauss linking invariant. 
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I. INTRODUCTION 

Subject of this work is the dynamics of a chain obtained by taking the continuous limit 
of a freely jointed chain consisting of N — 1 segments of length a and N beads of mass m 
attached at the joints between two consecutive segments. This problem has been addressed 
in the seminal paper of [jj using an approach based on the Langevin equation [31]. It was 
shown in |l| that the condition of fixed length segments becomes in the continuous limit a 
constraint which is similar to that of incompressible fluids in hydrodynamics. The authors of 
Ref. [l| have also described several interesting regimes in which their model of a constrained 
chain can be applied, like for instance an isolated cold chain or a hot polymer in the vapor 
phase. The statistical mechanics of a freely jointed chain in the continuous limit has been 



later investigated exploiting different methods, see for example 



Interesting related 



results may be found also in Ref. [6| . Up to recent times, however, most of the developments 
in the dynamics of a chajnwith rigid constraints have been confined to numerical simulations, 
see for example Refs. 



3,y,y. 

To overcome at least in part the complications of the dynamical case, it has been proposed 
in Ref. [l0( a path integral framework for the dynamics of the constrained chain discussed 
in l|. The resulting model, which describes the fluctuations of a two dimensional chain, 
is a generalization of the 0(2) nonlinear sigma model. For this reason, it has been called 
generalized nonlinear sigma model or simply GNLcrM. The relation of the GNLcrM with 
the Rouse model H] has been studied in details in Ref. [id ]. A difference between the 



two models concerns the scales of time and length at which the chain is observed. In the 



Rouse model only the long time-scale behavior of the chain is considered [12j |. On the other 
side, the GNLoM takes into account the short time-scale behavior and the finest details 
of the chain. These facts make the GNLcrM suitable to study the response of a chain to 
mechanical stresses in micromanipulations, for instance when it is stretched under a constant 
force |l3l . Il4j . Indeed, some experiments point out that the freely jointed chain model is 
able to capture the behavior of DNA in the limit of low applied forces [15( . 

The GNLcM does not take into account the hydrodynamic and self-avoiding interactions. 
The lack of hydrodynamic interactions limits its validity to the cases in which the motion 
of the beads is slow. This happens for instance when the viscosity of the fluid is large or 
the temperature is low. The conformations of the chain change slowly also in the presence 



of stiffness. The treatment of chain stiffness, a feature which was missing in the formulation 
of the GNLcrM of Ref. [10J], will be included in this work. A concrete application of the 
GNLaM could be polymers in a very dilute solutions at the so-called point, in which self- 
avoiding interactions play no role. The assumption that the chain is phantom, i. e. it can 
cross itself, is however dangerous at the G point because in that case it is very likely that the 
chain is knotted [l6[] and one should take into account the resulting topological constraints. 
In general, the fixing of constraints in (stochastic) dynamics requires some mathematical 
effort [7|, ll2|, ll7|, ll8( . The field theoretical formulation of the chain dynamics provided by 
the GNLcrM has the advantage that it is relatively easy to add further constraints, like for 
instance those which are necessary to impose topological conditions in the case of ring-shaped 
chains. 

The main goal of this work is the development of possible applications of the GNLcrM 
model. The most important result is indeed the calculation of the expectation values of 
two observables in a semiclassical approximation, which is valid if the changes in the chain 
conformation due to the fluctuations are small. This may happen when the chain is relatively 
rigid or in the following two cases: The temperature is low or the chain is moving in a 
very viscous solution. All these situations are compatible with the conditions of validity 
of the GNLcM mentioned before. The first observable which we consider is the dynamical 



form factor of the chain 12J. The second observable is a straightforward generalization 
to dynamics of the static form factor. It is shown that this observable is related to the 
average distance between two points of the chain. The calculation of both observables is 
complicated by the presence of ultraviolet divergences, which are regulated with the help of 



the zeta function regularization [19|. Let us note that divergences of this kind do not appear 
in analogous computations of the dynamical form factor performed using the Rouse model 

Another purpose of the present work is to improve the formulation of the GNLcrM given 
in [lO], making it more suitable for concrete applications. For this reason, we consider here 
the dynamics of a chain in d spatial dimensions. This case leads to a GNLcrM with 0(d) 
group of symmetry, which is a straightforward generalization of the two dimensional model 
already discussed in [10]. With respect to Ref. (lOj, we have also included in our approach the 
bending energy of the chain. In order to make the description of the chain dynamics closer 
to realistic situations, a method to take into account the topological entanglement of two 



closed chains is proposed. The topological constraints are imposed using the Gauss linking 
invariant. Unfortunately, it is not possible to apply to dynamics in a straightforward way 



the strategy based on Chern-Simons field 



of polymers, see for example Refs. 21 



22 



;heorv which is used in the statistical mechanics 
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261 ] . The main difference from statics 



is that in dynamics one has to take into account the motion in time of the chain. This 
implies that, rather than with the one dimensional trajectory of the chain, one has to 
deal with the two dimensional surface that the chain spans during its motion. To cope 
with this situation, we have generalized the multi-component Chern-Simons field theory 
of statistical mechanics to four dimensions. Mathematically, it is not possible to do that 
while keeping the topological invariance of the theory with respect to diffeomorphisms which 
depend both on time and on the spatial dimensions. However, the condition of invariance 
under diffeomorphism depending on time is not strictly necessary in the case of a non- 
relativistic chain and has been relaxed. 

The presented results are organized as follows. In Section II the problem of the dynamics 
of a chain in d dimensions is mapped into an 0(d) GNLcrM. The generating functional of the 
correlation functions of the bond vectors is expressed in the path integral form. In Section 
III the background field method is applied to the computation of the generating functional. 
Particular care is dedicated to the boundary conditions imposed on the fields to allow the 
freedom of performing integrations by part in the action without producing unwanted and 
cumbersome boundary terms. The action of the GNLaM is modified in order to take into 
account the bending energy of the chain. In Section IV the dynamical form factor and 
another related observable are computed. In Section V a model of two entangled chain is 
presented. Finally, our Conclusions are drawn in Section VI. 

II. FORMULATION OF THE MODEL 

In this section a path integral formulation of the dynamics of a freely jointed chain of 
length L is provided. The chain is regarded as a set of N beads connected together by N — 1 
segments of fixed length a. In addition, N, L and a satisfy the relation L = Na. Denoting 
with R„(£), n = 1, . . . , N the positions of N beads, it is possible to describe the fluctuations 
of the chain as a random walk of the beads constrained by the conditions: 

IR^-R^x^l^a 2 n = 2,...,N (1) 



These conditions are required by the fact that the length of the N — 1 segments connecting 
the beads is equal to a. We also demand that at the initial and final instants t — and t = tf 
the n-th bead is located respectively at the positions R„(0) = Ro, n and R n (£/) = R/,n- At 
this point, following Ref. lO], we introduce the probability function ipN which measures the 
probability that the chain after a given time tf passes from an initial configuration Ro, n to a 
final configuration R/, ra . Using an approach which is widespread in the statistical mechanics 
of polymers subjected to constraints, we define ipN as follows: 

* = /.,„,,, BR.W ■ ■ ■ /»„,.„=»,,„ MM«) «p {- E /„" «^} 

R 1 (0)=R 0jl R JV (o)=R 0>n *. n=X ) 

X TT S f iRnffl-Rn-lWI 2 _ ^ (2) 

n=2 V « / 

where .D denotes the diffusion constant. 

The path integral (j2j) describes the random walks of the A^ beads composing the chain. 
The insertion of the Dirac delta functions is needed to enforce the conditions (1T1), which 
describe the rigid constraints due to the non extensibility of the individual segments. We 
remember that the diffusion constant D satisfies the relation D = fiksT, where \x is the 
mobility of a bead, ks is the Boltzmann constant and T is the temperature. Moreover, 
a = — , where m is the mass of the bead and r is the relaxation time which characterizes 
the ratio of the decay of the drift velocity of the beads. Supposing that the total mass of 
the chain is M, we have of course that m = M? = ^a. Thus, Eq. (E|) becomes: 

N 



t 



N 



n=l' R n (0)=R 

N 



" r { M N /■*/ 



R„(i)-R„-i(t)r 



*R f [ ^ - 1 ) (3) 

The limit A^ — > oo, a — > in which the continuous chain is rigorously recovered has been 
already discussed in Ref. [10| in the two dimensional case. The extension to d dimensions 
is straightforward. Basically, the continuous limit consists in the following replacements of 
the basic ingredients appearing in the path integral of Eq. ([3]): 

N 



II /.»cv)=h,,» m ^ ^ JvR(t,s) 

' ~*-n,(0)=B. ,n 

Va/' dtRl(t) — ► [ ' dt [ ds~R 2 (t, s) 
~, Jo Jo Jo 



g s f \R n (t) R n -!(t)| 2 _ \ _^ ^^ s) _ 1} (4) 



n=2 \ « 



R/ >n — > R/(s) 

R(),n ► Ro(s) 



where s is the arc-length of the chain and < s < L. We have also introduced the notation 
as 



R' = Ip. Applying Eqs. (jl]) to Eq. Q, the probability function ip N becomes: 



*(%(*), Ro(a)) = J R(tfs) ^ f{s) VR(t,s)V\(t,s)expl-cJ o f dtJ Q dsR 



R(0,s)=R (s) 

xexpli f f dt f ds\{R' 2 - 1) I (5) 

with c = Ak m TtL ■ In the above equation the Lagrange multiplier A = X(t, s) has been 
introduced in order to represent conveniently the functional Dirac delta function appearing 
in the right hand side of Eq. (jlj) . 

Formally, the path integral in the right hand side of Eq. (jSJ) resembles the partition 
function of a quantum mechanical chain with constant density mass ^ after the analytical 
continuation to purely imaginary times: 

S = ^r dt dsB 2 (6) 

2L jo Jo 

To stress the close analogy with quantum mechanics, we remark that in Eq. ((Sj) the action 

So is multiplied by the inverse of the factor k = 2ksTr. It is known that k plays in 

the Brownian motion the same role of the Planck constant, due to the well known duality 



between quantum mechanics and Brownian motion 27|. One may show that the action So 

originates from the continuous limit of the kinetic energy of a free chain, see Ref. 10J in the 

two dimensional case and Ref. [28] in three dimensions. As we see from Eq. (jHD, the presence 

of rigid constraints is responsible for the appearance besides the action So of an additional 

nonlinear term given by 

Si = -i f tf dt f ds\{R' 2 - 1) (7) 

Jo Jo 

The Lagrange multiplier X(s, t) in Si closely resembles the pressure in incompressible hydro- 
dynamics, as it has been noticed in [l|. It expresses the fact that the segments composing 
the chain have a fixed length and thus they may not be compressed. 

To conclude this introductory Section, we specify the set of boundary conditions satis- 
fied by the bond vector R(t, s) in the probability function (jSJ). First of all, the boundary 



conditions at the initial and final instants and tf are given by 

R(0,a)=Ro(s) R(t f ,s) = R f (s) (8) 

where Ro(s) and R/(s) are static chain conformations. Additionally, it will be convenient 
to choose boundary conditions with respect to the arc-length s which allow integrations by 
parts in this variable without generating cumbersome boundary terms in the actions So and 
Si. To this purpose, we will limit ourselves to the following two choices: 

1. periodic boundary conditions in the case of a ring-shaped chain 

R{t,s + L) = R(t,s) (9) 

2. fixed end boundary conditions in the case of an open chain 

R(t,0) = ri R(t,L) = r 2 (10) 

where ri and r 2 are constant vectors. For the Lagrange multiplier A one may impose trivial 
boundary conditions in time 

A(0,s)=A(t/,s) = (11) 

and boundary conditions analogous to Eq. (Q in the case of a ring-shaped chain. For an 
open chain it is sufficient to require that: 

A(t,0) = X(t,L) = (12) 

III. THE GENERATING FUNCTIONAL *[J] AND THE BENDING ENERGY 

To the probability function \I/(R/(s),Ro(s)) of Eq. (J5]) we associate the following gener- 
ating functional \I/[J]: 

*[J] = / VRV\e' 1 ^ So ' Sl e^ fdt Io dsJ - R (13) 

Jb.c. 

where Sq and Si have been defined in Eqs. (jSJ) and (JTj) and J is an external current. The sub- 
script b.c. near the integration symbol means that an appropriate set of boundary conditions 
among those of Eqs. (l8])- (fT2l should be chosen. If J = 0, one obtains back the probability 



8 

function \l/(R/(s), R (s)) of Eq. (J5J). In the right hand side of Eq. (fT3l) we perform the 
following shift of variables 

R(t,s) = R b (t,s) + R q (t,s) (14) 

The field R q (t, s) describes the fluctuations around a fixed background chain conformation 
R b (t, s). We require that R b (t, s) satisfies the boundary conditions 

R 6 (t /)S ) = R / (s) 
R b (0,s) = R (s) 

with respect to the time t. The boundary conditions in the variable s corresponding to 
Eqs. © and (HOD are: 



1. ring-shaped chain conformations 

R b {t,s + L) = R b (s) (16) 

2. open chain conformations 

R 6 (*,0) = r 1 R b {t,L) = r 2 (17) 

Finally, we demand that the background conformation fulfills the constraint 

R* = 1 (18) 

The fluctuations R 9 (t, s) obey instead the following boundary conditions with respect to the 
time variable t: 

R q (t f ,s) = R q (0,s)=0 (19) 

In the case of the variable s, the analogs of Eqs. ffT6|) and (IT7j) are respectively: 



1. ring-shaped chain conformations 

R q (t,s + L) = R q (t,s) (20) 

2. open chain conformations 

R,(£,0) = R q (t,L) = (21) 



The expression of \I/[J] in terms of R& and R g is 

^[J] = e -S 6 f VR ^ VXe -f*fdtf o L ds[cilt-i\{Il'l+2R q -Il> b )+Il q -(J-2cR b )] ^ 

Jb.c. 

In Eq. (|22p we have introduced the notation: 



dt I ds 
o jo 



cR£ + R 6 ■ J (23) 



Let us note that it is not necessary that the background field Rb satisfies the classical 
equations of motion related to the action ik So T + Si appearing in Eq. (TT3|) . This would be a 
severe restriction. We recall in fact that the compatibility with the constraint ffl8l) requires 
that the solutions of the classical equations of motion R^(s) are static chain conformations 
independent of time 10(] . The main disadvantage of static background conformations is that 
with this choice the initial and final conformations of the chain in the boundary conditions 
(TT5|) must be the same, i.e. R/(s) = Ro(s) = R^(s). On the other hand, we have seen here 
that more general background fields can be considered. Their only effect is the addition to 
the external current J of the term 2cRb, see Eq. ( 1221) . Of course, for static solutions R& = 
and this term vanishes identically. 

Let's now investigate how it is possible to include in our approach the stiffness of the chain. 
To this purpose, it is convenient to require that R& is a static background conformation of 
the chain. Under this hypothesis, Eq. ([T4l becomes R(t, s) = R&(s) + R g (t, s). Taking 
the derivative with respect to s of both members of this equation, one obtains H.'(t, s) = 
R' b (s) + R' (£, s). At this point, only the fluctuations are responsible for the time variation 
of the vector field R'(t, s) which is tangent to the chain trajectory for each value of the 
arc- length < s < L. In fact, the contribution of the background vanishes identically: 
R^s) = 0. Therefore, we may use the module of Rg(t, s) as a measure of how fluctuations 
are effective in bending the chain. Accordingly, we modify the generating functional ^[J] 
in Eq. (122]) adding a bending energy term as follows: 



^ Q [j] =e -s 6 f p RqPAe -/o /d '/o Ld 4 c ^ +5E ^ R '^ iA(R?+2R ^ R ' 6)+R ^ (J_2c ^ ) ] (24) 

Jb.c. 

The term added to the action of the fields R q and A is: 

It is easy to realize that the parameter a has the dimension of an energy per unit of length. 
Thus, H a {t) = a J dsH'g represents the bending energy measured in energy units of IzbT 
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which thermal fluctuations deliver in the unit of time r to the system. To have stiff chains, 
the fluctuations of the tangent vectors must be small, i.e. 

R; 2 < 1 (26) 

This situation is verified in the following cases: 

1. 2ksTr is small 

2. a is big 

The first condition implies that either the temperature is so low that fluctuations become 
negligible or that the chain is fluctuating in a very viscous environment. In both situations 
one expects that conformational changes due to thermal fluctuations are small, so that the 
chain may be considered as rigid. The interpretation of the second condition is straightfor- 
ward. Finally, we notice that it is not possible to substitute naively R' with R' in Eq. (|25|) 
because in this way the bending energy term would be trivial due to the rigid constraints 
which require that R' 2 = 1. 

IV. COMPUTATION OF THE DYNAMICAL FORM FACTOR AND RELATED 

QUANTITIES 

This Section is devoted to the computation of the average values of two physically in- 
teresting observables. The average (O) of an observable O will be evaluated using the 
distribution: 

Jvp(K q ,X) =JvR q Jv\e-fo fdt foH^-^ 2 + 2 ^M (27) 

The boundary conditions satisfied by the fields R g and A are those specified in the previous 
Section. Let us note that with this choice of boundary conditions T>p(Ti q , A) is exactly the 
measure appearing in the path integral of the generating functional of Eq. (1221) as it should 
be. In writing Eq. ( 1271) we have required that the background field R& satisfies the classical 
equations of motion, so that R& = 0. Without this condition, it is very difficult to perform 
analytical calculations, even in the stiff chain approximation of Eq. ( I26i) . We have also 
neglected the stiffness term S a of Eq. fT25|) and the irrelevant normalization factor e~ Sb . 



11 

First of all, we consider the quantity: 

where 

&(t, s) = zk [S(t - t 2 )5(s - l 2 ) - 8{t - t!)S(s - h)] (29) 



In the above equation k is a constant vector and we have assumed that < t\ < £2 < 



tf- 



' [ L dh [ L dh [Wo' <* I? *****&)] = i r * /" ^2 (e lk -(%W-%W)\ ( 30 ) 
jo Jo L J L Jo Jo x ' 



Let us note that the observable (128j) is related to the dynamical form factor, see Ref. 12] 
for an introduction to that quantity. As a matter of fact, substituting the current ( 1291) in 
Eq. (j28p and integrating over l\ and Z2 we obtain: 

L 2 

where we have added the normalization factor j?. The quantity in the right hand side of 
the above equation is nothing but the dynamical form factor of the chain. 

We are now going to compute the expression of \l/(£i). Looking at the integration measure 
of Eq. (1271) . it is easy to realize that ^(£1) coincides with the generating functional *&[J] 
in the special case in which J is the external current £1 of Eq. f[2"9l . As we will see, the 
presence of Dirac delta functions in £j produces ultraviolet divergences in the expectation 
value \P(£i) which should be properly regulated. For simplicity, we will consider here stiff 
chains. As explained above, see Eq. (T26]) . this means that the changes due to the fluctuations 
of the vectors tangent to the chain's trajectory are relatively small, so that it is possible to 
neglect in the probability distribution ( 1271) the quadratic term in R' : 

r; 2 + 2r; • r; ~ 2r; • r; (31) 

Taking into account Eq. (l3~Tj) . the expression of ^(£1) may be approximated as follows: 

*(&) ~ [vR q Jv\e-Jo fdt fo ds l c ^ + y R <} (32) 

where 

j = £i-i^(AR' b ) (33) 



After a straightforward gaussian integration over R g in Eq. (1321) . we obtain 

#(£1) = fv\e sw (34) 
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with 



S(X) = - f tf dtdt' f ds)(t,s)G(t,t')tit',s) (35) 

4 jo jo 

In the above formula G(t, t') denotes the Green function satisfying the equation 

d 2 G(t,t') 



2c- 



dt 2 



-5(t,t') 



(36) 



The boundary conditions of G(t, if) at both initial and final instants t — and tf = are the 
same Dirichlet boundary conditions of the fields R a given in Eq. (TT9I) . The Green function 
G(t,t') may be written in closed form as follows 29]: 



G(t,t>) = -- 



*('-'/), (t _ + '(*-'/: 



9{t' - 1) 



(37) 



Here 9(t) is the theta function of Heaviside. Later on it will be necessary to evaluate G(t, t') 
on the line t = t'. In this case the right hand side of Eq. (1371) becomes proportional to 9(0), 
which is not a well defined quantity. This is rather a problem of the chosen representation 
than an intrinsic flaw of the solution of Eq. (l3Tj) . For this reason, it will be useful to derive 
a series representation for G(t, t'), which is regular when t = if . To this purpose, we use the 
definition of a Green function in terms of its eigenvalues and eigenfunctions: 

fn(t)fn(t') 



G(t,t') 



E 



A, 



where the eigenfunctions / n (t)'s satisfy the equation 



It is easy to show that 



and that 



2C QpUt) = Xnfn(t) 



m- $***=*■ 



A, 



2n 2 n 2 c 



(38) 

(39) 
(40) 
(41) 



where n > in both Eqs. ff40|) and (j4Tj) . Inserting Eqs. ff40l and ff4T|) back in fl38|) we obtain: 

tf 



g(*,0 = E 



en 2 ir 2 



nir(tf — t) mr(tf — t') 
sin sin - 



ti 



t, 



(42) 



n>0 °" " l / ''Z 

Remembering the definition of the current j of Eq. (1331) . it is easy to realize that the action 
S(X) of Eq. (j55|) is gaussian in the Lagrange multiplier A: 



S(\) = - f tf dtdt' f dsG{t, if) 
2 Jo jo 



Ci(t, s) • 6(^', s) + iX(t, s)K' b (s) ■ 1L&V, s) 



+ iX{t', s)R' b {s) ■ —^{t, s) + X(t, s)R' b (s) ■ g-^{X{t', s)R' b {s)) 



(43) 
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At this point we require that the background conformation Rft(s) describes a ring-shaped 
chain placed in a two-dimensional subspace. For instance, we may choose 



R&(s) = / du (cos (p(u), sin((/?(ti)), 0, ..., 0) + ~R b:0 
Jsq 



(44) 



So is the arc-length of the point R^o belonging to the background conformation. Clearly, 
the above expression of the background field R&(s) satisfies the constraint (Ti~8l) and the 
periodicity conditions ( 1T61) provided the function (p(u) is periodic modulo 2ir: <p(u + L) — 
<p(u) + 2kn, k — 0, ±1, ±2, . . .. Using the ansatz (14*41) it is easy to show that in the action 
S(X) of Eq. (143]) the last term in the right hand side may be rewritten as follows: 



1 /•*/ 

2 Jo 

1 /■*/ 



5 2 



dtdt' f dsG(t,t')X(t,s)R' b (s) ■ —(\(t',s)R' b (s)) - 
Jo os z 

dtdt' f dsG{t, t') [\{t, s) (\"(t', s) - \{t', s)<f' 2 {s) 



2. Jo Jo L v /J 

In deriving Eq. (1451) it has been exploited the fact that R' b is an orthonormal vector satisfying 
the relations R' b = 1 and H' b ■ R' b = 0. Thanks to Eq. ( 1451) . the path integral over A ( |34l) 
becomes: 

#(£0 = /p^ e |/ t/ ^t'/ Z 'd s G(t,t')[^(t lS ).Ci(t', S )+iA(t, s )^,T(t', S )+iA(t', S )6,T(t,s)] 

x e -5/o t/d ' d4 7o LdsG ( t ' i ')[ A '( t ' s ) A '( i '' s )+^ 2 W A (*' s ) A ( t '' s )] (46) 



where we have set 

Z 1>T (t,s) = R' b {s)-g 1 {t,s) 

After a straightforward integration over A, we obtain: 



*(& 



3 Ji+J 2 



with 



and 



Ji = - [ f dtdt' [ dsG(t,lf)$i(t,s)£i(lf,s) 

2 Jo Jo 



I 2 = -- f * dtdt! f dsds'G{t, t')K(s, s% T (t, s)£i T {t', s' 
2 Jo Jo 

In Eq. (150]) K(s, s') denotes the Green function satisfying the equation 

^-(cp'(s)) 2 ~\k(s,s') = -5(s-s') 



(47) 



(48) 



(49) 



(50) 



(51) 
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At this point, we may proceed our calculation of ^(£1) considering general background 
conformations of the form (J441) . However, one should keep in mind that, if the function 
ip is too complicated, the solution of Eq. ( 15T1) is not known explicitly. For this reason, we 



will concentrate here on the particular case in which ip(s) 

conformation R&(s) has the shape of a circle of length L: 

_ , . L ( 2ns 2ns 

R H S ) = 2~ I cos ~F~' sm ~r~' °' " " " > ° 



2-ks 
L ' 



so that the background 



(52) 



Let's now evaluate the integrals appearing in the two exponents in the right hand side of 
Eq. ( j48l) . The explicit calculation of I\ and I2 will be performed in Appendix A. Only the 
final results are provided here: 

h = (53) 

and 



E 

a=l 



hi K 



L 



fhj Ki n 



f \ x — ^ i i / /7\ / /7\ "^ / \ x — ^ i i fi/i \ n / 1 \ 



- E { frk 



2 -- J 



»J— 1 



G(ti,t 2 



4vr 2 - 



5 2 K(/i,/ 2 ) , 



azidz 2 



4(y^,,(y + 






dh 



x bAh)xlAh 



+ dK{ h M) xlAh)x' b ^l 2 ) + K(l u h)xl{liKAh) + 



dU 



In the above equation we have put 




tf 



g(t f ,t a ,c) = Y, -^-^sin 



nn(tf — t a ) 



n>0 



cn^n' 



t. 



Lr\V a , T c 



(54) 



(55) 



K(s,s') = -8(s,s') 



(56) 



with a = 1, 2 and a being a constant defined in Eq. (1A19J) . 

The function K(s,s') appearing in Eq. (1341) is the Green function of Eq. (|5"T1) . If the 

background conformation is given by Eq. f[5"2l) . K(s,s') satisfies the relation: 

'(?_ 4vr2" 
ds 2 L 2 

An explicit expression of the solution of Eq. (|5"61) in the form of a Fouries series is given in 

the Appendix, Eqs. (1A13J) and (1A14D . Finally, in Eq. (1541) the components of the background 

conformation field H b have been denoted with the symbols x' bi (s), i = 1,2. The particular 

choice of R&(s) made in Eq. (1521) implies 

4,i( s ) = cos^(s) xl 2 (s) = sinv?(s) (57) 
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with if(s) = 2m. Let us note that the delta functions present in the external current 
£i(t, s) of Eq. (129]) are responsible for the self- interactions of the two points on the chain 
corresponding to the values of the arc-length s — l\ and s — l 2 - These self- interactions 
introduce infinities in both integrals I\ and I 2 . Such infinities have been regulated in order 

;o obtain the final result of Eqs. (1531) and (1541) with the help of a (^-function regularization 

19j. At the end it is possible to write: 

*(&) = e h (58) 

where I 2 is given in Eq. ( |54l) . 

In a way which is analogous to that used to calculate the quantity ( 1281) one may compute 
also the following observable: 

#(£ 2 ) = L- J* f dt f o L ds&(t,syR q (t,s)\ ^ 

where 

&(*, ^) = T^V^ - *i)0(*a - t) [S(s - h) - 6{s - h)} (60) 

12 — 1\ 

^(£2) provides a measure of the average distance between two points of the chain over the 
time t 2 — t\. However, ^(£2) ma Y also be used in order to estimate the distance between 
two points at any given instant t\. As a matter of fact, substituting the expression of the 
current (1601) in Eq. (159]) and taking the limit £2 — ► t\, it turns out that: 



e -/oWc>^)- R ^%(£ 2 ) 



e ik.(R(ti,J a )-R(ti,h))\ ( gl ) 



where Xi(t,s) denotes the i— th component of the vector R(t, s). Expanding the exponent 
in the right hand side of the above equation up to the second order we obtain 

*(6)| tl=t2 ~ (l+ik-(R{t u h)-'R.(tuh))- Y.hikjixiituW-Xiituh)) 

x(x j (t 1 ,l 2 )-x j (t 1 ,l 1 )) + ...) (62) 



It is easy to show that, for instance: 

d 2 



ak 2 



*(6 



*1=*2 
k=0 



R{t h l 2 )-R(ti,h)\ 2 ) (63) 



confirming the close relation of the observable ^ (£2) with the average distance of two points 
of the chain. 
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The computation of ^(£2) can be performed in a way that is analogous to the calculation 
of ^(^i). We report here only the result 



(*(&)> = exp 



2 
/ , hkjBij 



(64) 



where 



and 



A d* KM) A dK(l lt l 2 ) l \ 

+ A -r^-x' h ,i{h)xlj{h) - ^- (xl&iWbjih) + x'l^l 2 )x'lj{l 2 )) 

+7 (K(hj2)a%,i(h)x'lj(l2) + K(l 2 , JiK'^K'^Ji)) (65) 



/■*2 dt /-*2 dt' „, .. ,„„. 

A= / -/ -G(M') (66) 

Jti l 2 — X\ Jt\ i 2 — t\ 



V. A MODEL OF THE DYNAMICS OF TWO TOPOLOGICALLY 

ENTANGLED CHAINS 

In this Section we discuss the physically relevant case in which d — 3. The single chain 
model will be extended to two chains including topological interactions, which in three space 
dimensions become relevant, in particular when the chains are near the ©—condition. Let 
us consider two closed chains C\ and C 2 of lengths L\ and L 2 respectively. The trajectories 
of the two chains are described by the radius vectors Ri(t, Si) and TL 2 (t, s 2 ), where < s\ < 
L\ and < s 2 < L 2 . The simplest way to impose topological constraints on two closed 
trajectories is to use the Gauss linking number \\ 

X (t, C 1} C 2 ) = ^-i dR 1 4 dR 2 x ^ ~* 2 l (67) 

An Jci Jc 2 |Ri — K2I 

If the trajectories of the chains were impenetrable, then \ would not depend on time, since it 
is not possible to change the topological configuration of a system of knots if their trajectories 
are not allowed to cross themselves. However, since we are not going to introduce repulsive 
interactions between the two chains which could prevent their crossing, we just require that, 
during the time tf, the average value of the Gauss linking number is an arbitrary constant 
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m = T- I"' xfaCuCjdt (68) 

t f Jo 

Our starting point is the probability function of two free chains: 

¥(<?!, C 2 )=ff[ [VR t VX t ] e-^ 1)+s( ^ (69) 

where, in agreement with Eq. (jSJ), the actions S^ 1 ) and S' 1 - 2 -' are given by: 

S-W = / 4/ (it / t/ Si [cR 2 - zA<(Rf - 1)] i = 1, 2 (70) 

jo jo 

In order to add the topological interactions, we introduce in the above functional a Dirac 
^-function which imposes the constraint (1681) . In this way we obtain the new probability 
function: 



<MCi,C 2 ) = fY[[VR t VX t ]5[m- I* — i dRi-/ dR 2 x ^ "« J 

J jt=i \ Jo Aixtf Jci Jc 2 |Ri-R 2 | d / 

Xe -(5W+5«) (71) 



Exploiting the Fourier representation of the Dirac 5-function, the probability function 
^m(Ci,C 2 ) takes the form 

* m (C 1} C 2 ) = / dA^ A (d, C 2 )e-^ A (72) 



where 



MCi.Ca) = /n^RiPAje-^+^^^^^^^-^^^iS^ (73) 

•^ »=i 

At this point, after introduction the three dimensional spatial indices fi,u,p = 1,2,3, we 
state the identity: 

a f tf dt / m f m (Ri - R 2 ) 
A / & dRi • & dR 2 x y 



47ft/ JCx Jc 2 |Ri — R2 

+00 



<fy / d 3 x / dt/ / d 3 y Jf (77, xjG^fa, 77'; x, y) J 2 V, y) (74) 



In the above equation we have defined the following currents: 



Ml,*) = 7* f tf fKv-t) f Lt ds l ^p^5^^-R l (t,s l )) (75) 

jo £/ Jo asj 
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with 7i = 57- and 72 = A. G^rj, rf\ x, y) is the propagator of the field theory 

Scs = ~ d V d 3 xA^( V , x) • (V x x A® fa x)) (76) 

tf J— 00 J 

V x being the gradient with respect to the spatial variable x. Moreover, the A^(?7,x)'s, 
i = 1,2, are two vector fields defined in the Euclidean four dimensional space (77, x) and 
having three spatial components A®. Explicitly, G^^rj, rf\ x, y) is given by 



t f (x - y) 



p 



G^fa rf; x, y) = -J-e^p f- yJ 5{ V - 77') (77) 

Z7T |x — y|° 

Apparently, Scs is similar to the multi-component Chern-Simons field theory used to impose 



the topological constraints in the case of a static chain 30]. It differs however from it by 
the addition of the fourth dimension spanned by the coordinate 77, with —00 < r] < +00. 
This new coordinate is necessary to deal with the time variable t appearing in the dynamical 
case. Using the identity (1711) . it is possible to formulate the probability function \1>a(Ci ) C 2 ) 
of Eq. (1731) as a Chern-Simons field theory 

^ A (d, C 2 )= ff[ [VR t VX t VA^\ e -^cs+sw+s^) e ^ElJ^^J d ^Jt(rj,^\^) (7g) 

where the actions Scs, S^ 1 ', S^ have been defined respectively defined respectively in 
Eqs. ( 176]) and ( 17D1) . while the currents Jf(?7,x) are given in Eq. ( 1751) . 

VI. CONCLUSIONS 

The main goal of this paper was to make the GNLaM of Ref. [lOj more suitable to 
describe realistic systems and to compute the average of concrete physical quantities. In the 
introductory Section II, the continuous limit of a freely jointed chain consisting of A^ — 1 
segments of fixed length a and Af beads of mass m attached at the joints has been discussed 
in d— dimensions. The final model describing the dynamics of the continuous chain is a 
generalized nonlinear sigma model. The difference from the two dimensional case discussed in 
Ref. [lOj is that the underlying symmetry group is O(d) and not 0(2). This slight difference 
is enough to complicate the calculation of the generating functional y£f[J] of Eq. ( fl3l) even 
in the approximation of Eq. (f26l) . To obtain analytical results, one is forced to assume that 
the background conformations R;, are lying on a plane as it has been done in Eq. 
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In Section III we have introduced in our approach the notion of chain stiffness. The 
bending energy term S a of Eq. fl25l) has been added to the action of the GNLctM in Eq. ff24l) . 
Let us note that, due to the constraint R^ 2 + 2R' b • R^ = 0, S a may be treated as a linear 
term, in which the fluctuation R q is coupled to an external current proportional to R^'. 

The expectation values of two important observables have been derived in Section IV. 
The first observable ^(^i) is the dynamical form factor. The second observable ^(£2) ma Y 
be related both to the average distance between two arbitrary points of the chain at a given 
time t\ or to the average of that distance over a finite time interval. We have seen that the 
calculation of these observables is complicated by ultraviolet divergences, which have been 
cured using the zeta function regularization. The closed form of ^(^i) and ^(£2) in the 
approximation ( |26l) has been presented in Eqs. ( J54|) -(l58l) and (|6l j) -f l66l) respectively. 

Finally, we have proposed in Section V a way to describe the topological relations between 
two ring-shaped chains via the Gauss linking number. This result generalizes to dynamics 



the treatment of topological constraints presented in the case of statistical mechanics in 30]. 
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APPENDIX A: PROOF OF EQS. ([53]) AND (l54h 

First we evaluate the integral ii of Eq. (H9l) . which we rewrite here for convenience: 

I 1 = \ f tf dtdt' I dsG(M')6(M)£i(0) (Al) 

2 Jo Jo 

where £i(£,s) is given in Eq. (|29|) . The only potentially non-zero contributions come from 
the self- interactions of the two points located at arc- lengths s = l\ and s' = l 2 : 

h = -— J2 I' 1 dtdt' [ L dsds'G(t, t')5(t - t a )S(t' - t a )5(s - s')5(s - l a )8(s' - Q (A2) 
2 ^[Jo Jo 
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There are in principle other two contributions which are proportional to 8(li — I2) and thus 
vanish identically, since l\ 7^ li- The time integrations in Eq. (1A2|) do not pose particular 
problems. Using the prescription (1420 to evaluate the Green function G(t, if) at coinciding 
points, we obtain 

h = ~ E 9(tf, t a , c) [ L dsds'8(s - s')S(s - l a )8(s' - l a ) (A3) 

Z (3=1 J0 

Here g(tf,t a ,c) is the series given in Eq. (|55[ ). Unfortunately, the integrals over the arc- 
length s and s' are divergent and require regularization. To this purpose, we first expand 
the periodic ^-function 5(s) in Fourier series : 

+00 2ttikJ- 

S{*)= E ~ ir (A4) 



K = — OO 



After some calculations, it is possible to show in this way that: 

r L 1 +°° 

/ dsds'5(s - s')5(s - l a )6(s' -l a ) = T Y 1 (A5) 

7 ° L K =-oo 

Of course, J2t=-oo 1 ^ s divergent if left without treatment. To remove the singularities, we 

will use the zeta function regularization. This kind of regularization is based on the Riemann 

C-f unction: 

+00 -I 

COO = E "I (A6) 

and on the fact that, in the sense of the analytic continuation, one may write the following 

formal identity: 

+00 

C(0) = E 1 (A7) 

K=0 

On the other side, Eq. (1A7J) implies that: 

+00 

E 1 = 2C(0)-1 (A8) 

K= — OO 

Applying Eqs. flA7[) and (|A8[) to Eq. (|A5[) and substituting the result in the expression 
of I\ given in (IA30 . it is easy to realize that the integral /1 becomes: 

/i = -^/,t a ,c)E 7 (2C(0)-l) (A9) 

2 a=l L 

After an analytic continuation of the function Q(s) to the point s — 0, one finds that £(0) = |. 
Substituting this value of £(0) in Eq. (1A9I) . we obtain 

/1 = (A10) 
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This completes the proof of Eq. (I53l) . As expected, the self- interactions of the points at 
s — l\ and s — I2 with themselves do not give any contribution to ^(£1)- 

The situation is more complicated in the case of the second integral 1% of Eq. (1501) : 

J 2 = -I f h dtdt' f dsds'G(t, t')K(s, s')& r(t, s)6 T (t', s') (All) 

2 Jo Jo 

Using the definition (T4"T|) of the current £1^ one obtains: 



2 2 ^fc. L 

h = J2 JL-^^tftt^c) I dsds' 







K^s'W^x'^s') + ^f^<^K;(*0 



dK(s,s') , „ , dK(s,s') „ , , 

-HM x b,A s ) + — ^ — x bA s ) x b,j( s ) 



(9.S 



ds' 



8(s - l a )S(s' - l a ) 



L.7=l 2 ^ 



5 2 i^(s S 1 ) 



, 9K(s, s') , „ , dK(s,s') „ , , 

<9s X b,i\ S ) X b,j\ S 1 ' a 1 x b,i\ s ) x b,j\ s ) 



dsds' 
8(s - l a )5(s' - Q 



+ 



ti <-> U 



(A12) 



/1 <-> ^2 

To write Eq. (IA12[) . we have made some integrations by parts in the variables s and s'. These 
are allowed because of the choice of the boundary condition and of the fact that the current 
£1 vanishes at the boundary: £i(£, 0) = £(t,L) = 0. It is not difficult to show that, for 



symmetry reasons, 
of K(s,s'): 

where 
we obtain: 



dK(s,s') 
ds 



s — l a 

,' = lr 



0. As a matter of fact, using the Fourier representation 



K(s,s')= J2 e^ s -^K{n) 



K{k) 



dK(s,s') 



ds 



4vr 2 k 2 + 1 
1 +00 



ft- 



Analogously, one may show that — » , ' 



2vr w f- k> + 1 



0. 



(A13) 
(A14) 
(A15) 



There is only one term in the expression of I2 which is divergent and needs regularization. 
This is given by: 



2 2 



K; j h -j 



i J=l a=l 



L ^j*KM^ 



h, S ing =EE -hr9{tf,tonc) / dsds' ' — x' bi (s)x' bj (s')5(s - l a )5(s' - l a ) (A16) 



dsds' 
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Exploiting the Fourier representation (1A13J) of K(s,s'), it turns out that I 2)S i ng m ay be 
rewritten as follows: 



'2, sing 



2 2 hk- 

J2 J2^r9(tf,t a ,c)x' b!i (i a )x' bd (i a ) T J2 

i,j=l a=l 



1 +2? K 2 



L K= _oo K z 



(A17) 



-OO K 



T,t=-oo (l - m) and the fact that E+T-oc 1 = 



Applying also the identity Efcl-oo ^rpi 

as we have previously seen, we arrive at the final result in which the singularity of lining has 

been regulated: 



a 



l2,sing — ~J /^ L^i ~^~9{tf i^a, C)X bi {l a )X b j{l a ) 
^i,j=la=l A 



In the above equation we have put: 



(A18) 



a 



E 



(A19) 



ho= J:J:^g(t f ,t ai c) f dsds>xl(s)xl(s')5(s-l a )6(s' -l c 

i,j=la=l 



(A20) 



Let's now simplify the following term contained in I 2 : 

2 2 /••/■ ,-L 

After the integrations over s and s', one obtains from I 20 an expression which is proportional 
to the Green function K(s, s') computed at coinciding points s = s' = l a . Exploiting the 
Fourier representation (1A13I) . it is possible to check that K(l a , l a ) is convergent and is equal 
to: 



J\ yla, t a ) 



4^ 



where a is the constant given in Eq. ( 1A19I) . Thus, we may write: 

2 2 kk- nJ 

J 2,o = 2^ 1^ ^-9{tf,t a ,c)—x b!i (l a )x b!J {l a ) 

ij=l a=l / 47T 



(A21) 



(A22) 



All the other terms present in I 2 are divergenceless. At the end, remembering the expressions 
of the contributions I 2)S ing and J 2 ,o to I 2 given in Eqs. (1A18I) and (1A22I) respectively, we obtain 
the final result: 



2 ££. 



^ (M2) <^K,(* 2 ) + ^§^<^K,(y 



dhdl 2 



dh 



+ 9Ki y 2) x'UhKjih) + K(l u l 2 )xl(h)xl(l 2 ) + 



m 2 



K h^ i 2 



(A23) 
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The above equation coincides exactly with Eq. f|54|) . 
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